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ABSTRACT. A semi-empirical m^ethod is used to characterize the 3s^3p^ — SsSp'' J=2 transition array 
in P II. In this method. Slater, spin-orbit, and radial parameters are fitted to experimental energy lev- 
els in order to obtain a description of the array in terms of LS-coupIing basis vectors. The various 
IC and CI amplitudes resulting from this model are then used to predict the branching fractions of 
transitions within the array. Results close to LS-coupling values are presented, and these are com- 
pared to branching ratios measured using beam-foil spectroscopy at the THIA laboratory. The work 
provides support for the hypothesis of Dr. Curtis that transition arrays with little upper state IC but 
significant upper state CI in atoms of low Z exhibit branching fractions close to LS-coupIed values, 
although the data are inconclusive in this respect. 
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1. Introduction 

In light atoms carrying a small central charge, 
relativistic effects upon the electron cloud are 
generally small. For this reason, it is con- 
venient to refer to levels within such sys- 
tems using the LS coupling model. However, 
in any multi-electron atomic system, LS cou- 
pling fails as a precise model due to direct 
and exchange coulomb interactions between 
various electrons, relativistic spin-orbit effects, 
and, somewhat less significantly, the magnetic 
dipole-dipole interaction between each electron 
and the nucleus. In order to take these inter- 
actions into account, the intermediate coupling 
(IC) model must be used. In this approximation, 
each electron is described by a distinct wave- 
function that is a linear combination of LS basis 
states with different L or 5, but the same J. The 
IC model preserves the convenience of utilizing 
LS selection rules w^hen identifying transitions, 
and it can be generalized to include configu- 
ration interaction (CI), in which the orthonor- 
mal basis set is extended to comprise LS states 
from different configurations. Because of this, 
the IC/CI model provides a useful picture of the 



atom that can be used to make predictions and 
isolate trends along isoelectronic sequences. 

In order to successfully use the IC/CI model 
to characterize an atomic system, it is nec- 
essary to obtain the mixing amplitudes that 
give the electronic eigenstates of the atom in 
terms of an LS basis set. These amplitudes 
may be determined ab initio, using iterative 
energy-minimization algorithms such as Multi- 
Configuration Hartree Fock (MCHF), or they 
can be arrived at semi-empirically ||9l. In this 
paper, we use the latter approach, which relies 
upon the high accuracy and precision of NIST 
energy-level data to obtain accurate results. 

2. The Semi-Empirical Method 

In the semi-empirical technique we use 
to model the system, the energies contained 
within the various interactions described above 
are treated as parameters and fitted to measured 
energy-level data. These parameters include the 
Slater direct and exchange energies F^ and G^, 
the spin-orbit energies Cfc/ arid the radial config- 
uration interaction energy R^. Linear combina- 
tions of the parameters are placed in a symmet- 
ric nondiagonal interaction matrix expressed in 
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the finite LS basis of interest, where each ma- 
trix element corresponds to the Hamiltonian be- 
tween two LS states. IC and CI mixing ampli- 
tudes are obtained by diagonalizing the matrix 
and finding its eigenvectors, which correspond 
to single-electron wavefunctions in terms of LS 
basis states. 

In order to diagonalize the matrix, w^e obtain 
a characteristic equation of the form 
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where each a„ is expressed in terms of the var- 
ious energy parameters. The roots of this equa- 
tion are the eigenvalues of the interaction ma- 
trix, w^hich would ideally correspond closely to 
the measured energy levels of the system. How- 
ever, due to the imperfect nature of the model, 
it is generally impossible to find real parameter 
values which give this result. Instead, w^e use 
the method of Z. B. Rudzikas [3j, noting that be- 
cause the characteristic equation can be factored 
into its experimentally measured roots, 

{X-Ei){X-E2){...){X-En) = 

the coefficients of A can be expressed as sums 
and products of experimental energies. For in- 
stance: [7] 
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Thus, the energy parameters may be adjusted 
so that the a„ in the characteristic equation fit 
these experimental coefficients. This is accom- 
plished using a numerical least-squares algo- 
rithm that searches parameter space for the 
global minimum of a x^ surface. When such 
a minimum is encountered, its parameter coor- 
dinates are printed and the interaction matrix 



may be expressed numerically. The eigenval- 
ues and eigenvectors of the interaction matrix 
can then be found, giving the IC and CI mixing 
amplitudes that facilitate semi-empirical com- 
putation of branching fractions and lifetimes. 



3. Application of the Semi-Emprical 
Method to J=2 Transitions in P II 

Here, we have applied this method to 
3s23p2 - 3s3p3 J=2 transitions in P II. Due to 
significant configuration interaction between 
the 3s3p^ and 3s^3p3d upper states, it was nec- 
essary to devise an 8x8 interaction matrix for 
J=2 that would take the levels of both of these 
upper configurations into account (see Figures 
1-3). This matrix included the average config- 
uration energies Ea and Eb in addition to the 
other parameters mentioned earlier, making ten 
parameters total. The coefficients of the var- 
ious parameters were determined from tables 
in Condon and Shortley IQI. The model was 
underdetermined, so in order make the system 
overdetermined and "steer" the x^ fit into a 
physical result with some empirical and ab ini- 
tio information, Ea and Eb w^ere set to the mea- 
sured average configuration energies, while Ri, 
Cp, and Cd were fixed to MCHIjjvalues: 



Parameter 


Value (cm ^) 


Ea 


78, 700 


Eb 


96, 500 


Ri 


13,348 


Cv 


295 


Cd 






Due to fixing these five parameters, the ac- 
tual parameter space w^as M^. The x^ fit 'was 
accomplished using Mathematica; however, the 
software's nonlinear regression package w^ould 
only fit a single analytic function to a set of data 
points, w^hile in our case it was necessary to fit 
eight distinct expressions simultaneously. For 
this reason, it was necessary to write out a man- 
ual weighted computation of x^ and then use 
the internal function NMinimize to search the 
parameter space (see Figure 4). 



Note that the p shell is half-filled in the 3s3p^ configuration, so {Hso) = 0, causing C,^ to appear only as a second- 
order, off-diagonal parameter in the interaction matrix. 



1: I 3s3p3 5S°) 2: | 3s3p3 ^D") 3: | 3s3p3 3p°) 4: | 3s3p3 ^D") 

5: I 3s23p3d ^D") 6: | 3s23p3d ^F") 7: | 3s23p3d ^F°) 8: | 3s23p3d ^D") 

Figure 1 . The ordered LS basis set of the J=2 interaction matrix 
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Figure 2. The J=2 interaction matrix 
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Figure 3. Diagonal elements of the J=2 interaction matrix 

chiSquared = 0; 

For[i = 1, i < 9, i + +, chiSquared = chiSquared + ( (energyCoef f icients [ [i; 

- Part [lambdaCoef f icients, 9 - i] ) '2 ) / (energyCoef f icients [ [i] ] "2) ] ; 

Figure 4. Manual computation of x^ in Mathematicaj 



To weight ^ for uncertainty, the code divides each squared difference by the square of the corresponding numerical 
coefficient of A. This is based upon the assumption that the uncertainties in the measured energies are essentially 
proportional to their magnitude. 
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The parameter search ultimately returned 
X^ ~ 2.76893 X 10"^^ and gave the following 
values for each of the five Slater parameters: 



Parameter 


X' fit (cm-^) 


MCHF (cm-^) 


F^(3p, 3p) 


9,646 


48, 076 


i^'(3p, 3d) 


48, 202 


5,334 


Gi(3p,3d) 


6,643 


2,824 


Gi(3s, 3p) 


59,117 


66, 845 


G3(3p, 3d) 


-22,040 


1,984 



While both of the G^ parameters exhibit order- 
of-magnitude agreement with MCHF values, 
both of the F^ parameters and the G^ param- 
eter differ significantly from the MCHF calcu- 
lation. However, this comparison may not be 
as unfavorable as it first appears. It is worth 
noting that at the MCHF parameter coordinates, 
the height of the x^ surface is 0.118162. Thus, 
while this value is many orders of magnitude 
larger than the best-fit value of 2.76893 x lO^^^, 
the 8x8 interaction matrix model is a good fit to 
measured energy-level data even with MCHF 
parameters. This is an encouraging conclusion, 
because MCHF is considerably more sophisti- 
cated than a nonlinear regression technique in 
this context and has been more reliably applied 
to complex atomic systems |5|. Therefore, al- 
though the 8x8 matrix model returns results 
considerably different from MCHF, its charac- 
terization of the system is approximately physi- 
cal. 

It is difficult to construct a conceptual link 
betw^een the mechanisms of a parameter-search 
procedure and the physical reality of an atomic 
system; however, some tentative comments 
may be made on these results. First of all, it 
is not certain that the disagreement between 
the MCHF results and the x^ results is en- 
tirely due to imperfection in the 8x8 model; it 
is possible for the 8x8 model to capture effects 
that MCHF does not. The apparent "switch- 
ing" of the F^ parameters is of particular in- 
terest in this respect. Basic physical intuition 
supports the MCHF result that F2(3p, 3p) > 
F^(3p, 3d), because it is reasonable to expect 
that (i?3d) > (i?3p), reducing the p/d coulomb 
interaction relative to the p/p interaction. Even 

4 This was suggested by Dr. ElUs. 



SO, the MCHF algorithm utilizes a central-field 
approximation that can leave out the electron 
correlation effects inherent in the actual multi- 
electron w^avefunction of the system. Thus, the 
X^ characterization of these parameters could 
indicate the presence of electron correlation and 
markedly distinct angular expectation values 
among the electrons of the 3p suborbitajj How- 
ever, the results in this regard are highly incon- 
clusive and such an interpretation represents lit- 
tle more than conjecture at this time. 

A comparison of fitted and experimental en- 
ergy levels provides another test of the x^ fit 
quality. In order to produce this comparison, 
eigenvalues and their corresponding eigenvec- 
tors were computed from the numerical J=2 ma- 
trix. LS states and eigenvectors were then as- 
signed to each other in pairs according to the 
largest component of each eigenvector. 



Config. Level 


^Fitted (cm ') 


i^Exp. (cm-^) [8J 


3s3p^ ^S^ 


45,665 


45,697 


3s3p3 ^D^ 


67,231 


65,272 


3s3p3 3po 


70, 689 


76, 764 


3s23p3d ^D^ 


81,792 


77,710 


3s23p3d ^D^ 


91,145 


104, 102 


3s23p3d ^F^ 


96, 173 


87, 804 


3s3p3 iD^ 


109, 604 


112,607 


3s23p3d 3p° 


111,287 


103,629 



This table shows that the largest disagreement 
between fitted and experimental energy levels 
is approximately 12,957 cm^^ (3s23p3d ^D^); 
however, the largest disagreement within the 
configuration of interest, 3s3p^, is a respectable 
6,075 cm~^ (^P2)- Thus, despite some flaws, 
the x^ model succeeded at describing the P II 
system to a precision commensurate with its so- 
phistication. 



4. Failure of the Semi-Emprical Method 
FOR J=l Ground Transitions in P II 

It is worth mentioning that similar applica- 
tion of the semi-empirical method to the J=l 
case for 3s^3p^ — 3s3p'^ transitions was not a suc- 
cess. A 7x7 interaction matrix w^as generated for 



3s3p3 ^D^)' = 0.975026 | 3s3p3 ^D^) - 0.219777 | 3s23p3d ^D^) + 0.031831 | 3s23p3d ¥^) 
I 3s3p3 3P^)' = 0.991192 | 3s3p3 ^P^) + 0.109902 | 3s23p3d ^P^) + 0.0694869 | SsSp^ ^D^) 



,3 3t-iO\/ 



3s3p^ ^D^/ 



3s3p^ °S^ 



0.963482 I 3s3p'^ '^D^) - 0.257754 | 3s^3p3d -^D^) 



3 5co\/ 



3 5co\ 



-0.999931 I 3s3p^ °S^ 



Figure 5. Significant components of eigenvectors for 3s3p J=2 level; 



levels in the 3s3p^ and 3s^3p3d upper states, but 
it was impossible to find a sufficiently low point 
in the x^ surface. As a result, an appropriate 
orthonormal eigenvector basis that would diag- 
onalize the matrix to measured energies could 
not be obtained. At the moment, it is not clear 
whether this was due to mathematical failure of 
the model or due to the presence of a signifi- 
cant interaction in the system that was left out 
of the calculation. The failure of the J=l model 
is unfortunate, and it is important to qualify the 
Slater parameters of the apparently successful 
J=2 model by noting that it is impossible to com- 
pare them against J=l values thus far. 

5. Computation of Semi-Empirical 
Branching Fractions for J=2 
Transitions 

As mentioned above, eigenvectors were ob- 
tained for the eight levels of the J=2 case af- 
ter the x^ fit produced Slater parameters. The 
eigenvectors described each IC/CI mixed state 
as a linear combination of LS basis states. Al- 
though the eigenvectors were found to have 
nonzero components along every element of the 
LS basis, each eigenvector had no more than 



configuration were written in terms of semi- 
empirical singlet-triplet mixing angles: 



'Pi) 

^P2) 



I 'Pi) 

cos ^2 



P2)-sin02 
sin 6*2 I ^P2)+cos6l2 



^D2) 
^D2 



In this case, 02 was found to be 1.630° [2J, and 
the J=0 levels of the low^er configuration were 
ignored because of the following El transition 
selection rules: IH 

A J = 0, ±1 (no to 0) 
AL = 0,±l(noOtoO) 

AS = 

The various line strengths from each upper state 
to each lower state were then computed (see 
Figure 8). Each line strength S^i was computed 
as the square of the reduced matrix element D^^l■. 
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where V is the dipole transition operator, Cl's 
four significant components. The most signifi- refers to an accessible componenlMof the primed 
cant components of each SsSp'^ eigenvector are lower state bra, A^s refers to a 3s3p'^ compo- 



displayed in Figure 5. The fit determined the 
3s3p^ ^S2 level to be essentially pure, so this 
level was ignored in computation of branching 
fractions. 

In order to compute branching fractions, it 
w^as also necessary to characterize the ground 
state appropriately, which was accomplished 
using a method formulated by Dr. Curtis jlj. 
In this method, the eigenvectors of the 3s^3p^ 



nent of the unprimed upper state ket \i^°j), and 
Bls refers to a 3s^3p3d component of the up- 
per state ket. The M factors are the LS-coupling 
transition elements, displayed in Figure 7. 

Each line strength was then corrected for 
wavelength, using: 
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All eigenvector components were used in the computation of branching fractions. 
That is, accessible by El transition rules to the upper state component in question. 



where the prime on ^^j indicates that this quan- 
tity is proportional, but not equal to, the transi- 
tion rate A^i- Given the corrected line strengths, 
it was straightforw^ard to compute the branch- 
ing fractions for each set of transitions: 



BF 
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ul 



ul 
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where the denominator on the right hand side 
is the sum of line strengths over all lower states 
that are El-accessible to the upper state u. 

Previously-measured experimental values 
for the branching fractions in question were un- 
available; however, it was possible to compare 



our semi-empirical branching fractions with the 
theoretical computations of S. S. Tayal and A. 
Hibbert. The 3p^ branching fractions were also 
measured at the THIA laboratory (see Figure 6); 
unfortunately, it was impossible to measure the 
3D2 branching fractions using THIA due to an 
especially large r, which is largely caused by 
configuration-interaction transition rate cancel- 
ing between the sp3 and spd components of that 
level. The various branching fractions and their 
comparisons are presented in the table below. 
Note: Tayal and Hibbert quoted /length values; 
these branching fractions were computed from 
those. 



Transition J Tayal BF [ilOJ Hibbert BF |6| LS BF Semi-Emprical BF Experimental BF 



0.233 ±0.0223 
0.767 ±0.0569 



3p^-3DO 
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3p,-3po 
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0.242 
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0.758 


0.748 


0.682 
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Figure 6. Measurement of 3s3p'^ 3po branching fractions at THIA. The intensities 
of the 3s23p2 3p^ ^^ 3p^ jj^gg ^^ 1305.48 A and 1310.70 A were measured using a 
Gaussian fit upon the sum of multiple spectroscopic scans. 
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Figure 7. The LS-coupling El transition elements, derived from Cowan by Dr. 
Ellis [9|. Note that 5*2 and the coefficient of fractional parentage cfp are given in 
the table at bottom. For computation of branching fractions, the hydrogenic ratio 
^pd = V'Vs'^sp was used, so that the radial integrals would cancel. 

For[i =1, i < 9, i++, 
For[j =1, j < 4, j++, 
For[k = 1, k < 9, k++, 
For[l =1, 1 < 4, 1++, 
If[l == KroneckerDelta[lowerSLJBasis [ [1, 1]], upperSLJBasis [ [k, 1]]] 
&& (upperSLJBasis [ [k, 2]] == lowerSLJBasis [ [1, 2]] || 

upperSLJBasis [ [k, 2]] - 1 == lowerSLJBasis [[ 1, 2]] || 
upperSLJBasis [ [k, 2]] + 1 == lowerSLJBasis [[ 1, 2]]) , 
relativeStrengths [ [i, j]] = relativeStrengths [ [ i, j]] + 

upperVecs [ [i, k] ] *lowerVecs [ [ j, 1] ] *mSelector [k, upperSLJBasis [ [k, 2]], 
upperSLJBasis [ [k, 1]], upperSLJBasis [ [k, 3]], lowerSLJBasis [ [1, 2]], 
lowerSLJBasis [ [1, 3]]]] 
]]]] 

Figure 8. Computation of all twenty-four 3s^3p^ — 3s3p^ J=2 line strengths in Mathematicaj 



The outermost loop cycles through all eight upper eigenvectors, the next loop cycles through all three lower eigen- 
vectors, and the inner loops cycle through the components of each upper-lower pair of eigenvectors. The "if" statement 
imposes selection rules and computes Dui- 
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6. Conclusion 

The theoretical, semi-empirical, and experi- 
mental values presented here are within 0.07 of 
pure LS-coupling branching fractions for these 
levels. This is to be expected, given the high 
(greater than 92%) leading components of all 
upper state eigenvectors in the x^ fit- This pro- 
vides support for Dr. Curtis' hypothesis that 
little upper state IC, but significant upper state 
CI, does not cause branching fractions to devi- 
ate heavily from LS coupling. 

For the branching fractions of the ^P2 up- 
per state, there is good agreement between the 
THIA values, Tayal, and Hibbert. Similarly, 
there is good agreement between Tayal and Hi- 
bbert for the ^D2 branching fractions. However, 
the semi-emprical branching fractions quoted 
here do not agree nearly as well with the other 



calculations. Due to the good agreement be- 
tween Tayal, Hibbert, and experiment, it seems 
likely that the semi-empirical approach is at 
fault here. Although the semi-empirical method 
produces branching fractions somewhat close to 
LS-coupling values, it appears to lack the so- 
phistication necessary to reproduce experiment, 
or the more advanced theoretical calculations 
for this system. 

The most interesting result of the semi- 
empirical method as applied to this set of tran- 
sitions in P II is the complete failure of the J=l 
model. Although this outcome could be due to 
limitations in the model, it could also indicate 
the presence of a heretofore unnoticed config- 
uration interaction or level crossing present in 
the system. Further study will be necessary to 
explore these possibilities. 
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